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Abstract 

We study the gluon propagator and the singlet potential in Landau gauge in 
the deconfmed phase of SU(2) lattice gauge theory, using both the standard Wilson 
action and a tree-level Symanzik improved action. From the long-distance behavior 
of correlation functions of temporal and spatial components of the gauge fields we 
extract electric (m e ) and magnetic (m m ) screening masses. For the magnetic mass 
we find m m (T) = 0.456(6) g 2 (T) T. The electric mass can be described by a next- 
to leading order ansatz, obtained from one loop resummed perturbation theory. 



However, the best description is given by m e {T) = y 1.70(2) g(T) T. The electric 
screening mass thus is different from its lowest order perturbative prediction even 
for temperatures as high as T ~ 10 4 T C . 



1 Introduction 



One of the main features of the high temperature plasma phase of QCD is the 
occurrence of chromo-electric and -magnetic screening masses (jn e and m m ) which 
control the infrared behavior of the theory 0. The electric screening mass, which 
has been known in leading order perturbation theory for some time (m e ~ gT), 
gives rise to the Debye screening of the heavy quark potential. The magnitude of 
m e influences strongly the existence or non-existence of hadronic bound states in 
the high temperature phase. An understanding of its temperature dependence is 
therefore essential for any further analysis of the quasi-particle excitation spectrum 
in the QCD plasma phase. Also, a non-vanishing magnetic mass, which is entirely 
of non-perturbative origin and generally is expected to be 0(g 2 T), does contribute 
in next-to-leading order to m e ||, [|. Beyond leading order perturbation theory one 
thus needs a non-perturbative analysis also for the quantitative determination of 
the electric mass^. 

The definition of gluonic screening masses is not without problems. In general 
the masses are related to the low momentum behavior of the gluon propagator, 
n w (po,p)- However, as the gluon propagator itself is a gauge dependent quantity it 
is not obvious that a "gluon mass" extracted from it would have a physical interpre- 
tation. Using the zero momentum limit (\p\ — ► 0) in the static sector (p = 0) of n MM 
does not yield a gauge invariant definition of the screening masses. It has, however, 
been shown || that the pole masses defined through m 2 ^ = 11^(0, \p\ 2 = —m 2 ^) 
are, within a wide class of gauges, gauge invariant to arbitrary order in perturbation 
theory. Exactly these pole masses are obtained from the exponential decay of fi- 
nite temperature gluon correlation functions at large spatial separations. Although 
these correlation functions have to be calculated in a fixed gauge the pole masses 
extracted from them will be gauge independent. Alternatively one may insist on 
defining gluon screening masses through gauge invariant operators, e.g. Polyakov 
loop correlation functions, which are related to the heavy quark potential at finite 
temperature. Observables, which project onto states with the correct quantum num- 
bers have been discussed in Ref. J?) and have recently been used to analyse the high 

a The non-perturbative structure of the electric mass beyond leading order perturbation theory 
has recently also been discussed in Q]. The interesting possibility has been raised there that 
non-perturbative effects may arise even without generation of a non- vanishing magnetic mass. 
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temperature phase of QCD within the framework of dimensional reduction || |9|] . In 
how far the definition of a screening mass through gauge invariant operators agrees 
with the pole mass has to be further analysed. The experience made with simi- 
lar concepts in gauge-Higgs models shows that there is the possibility that gauge 
invariant operators only project onto a superposition of several elementary gluon 
excitations iTDI . 



In this paper we analyse gauge dependent gluon correlation functions in Landau 
gauge. We have studied gluon correlation functions at finite temperature previously 
for the SU(2) gauge theory in a limited temperature interval (1 < T/T c < 20) Jll] 



This analysis will be extended here in various ways. We have calculated magnetic 
and electric screening masses at much higher temperatures, up to T~ 10 4 T c , in order 
to analyse the approach to the perturbative regime. Furthermore, we study this time 
the dependence of the gluon correlation functions on momenta. In addition, we now 
also analyse the color singlet heavy quark potential, which at high temperatures 
is dominated by single gluon exchange. Moreover, we have also used a tree-level 
Symanzik improved gauge action in addition to the standard Wilson action in order 
to get some control over systematic lattice discretization errors. 

In the next section we briefly discuss the observables we are going to analyse. 
In section |] we define the actions we have used in our simulations and describe 
the procedure followed to fix the relation between bare gauge couplings and the 
temperature. The determination of electric and magnetic screening masses from 
gluon correlation functions is discussed in section |j. In section || we discuss the 
determination of the electric screening mass from the color singlet heavy quark 
potential. Furthermore we describe the influence of the improvement of the action 
on this potential. Finally, we give our conclusions in section 10. 



2 Gluon Propagator on the Lattice 

The main purpose of this investigation is an analysis of the gluon propagator at 
small momenta. We will extract screening masses from correlation functions of gauge 
fields A^xq.x) at large spatial separations. We will analyse dispersion relations for 
these screening masses in the static (p = 0) sector. For this purpose we introduce 
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momentum dependent gauge fields, 



Mp ± ,x 3 ) = e ix ^A^x ,x ± ,x 3 ) , (2.1) 

xo,x ± 

and the corresponding correlation functions 

GJp±,x 3 ) = (ir A^p ± ,x 3 )AUp ± ,0)) (2.2) 



with x±_ = (xi,x 2 ) and p±_ = (pi,P2)- On a finite lattice, the momenta are given by 
Pi = 2irki/(aN i ), with ki = — ~iVj + 1, . . . , ~iVj and iVj being the length of the lattice 
in the i-th direction. 

The long-distance behavior of G yields the energy in the electric and magnetic 
sector respectively, i.e|] 

G e (p±,x 3 ) = G {p±,x 3 ) 

~ exp{-E e ( Pl _) x 3 } fora; 3 >l, 
G m (p±,x 3 ) = -(G 1 (p ±j x 3 ) + G 2 (p± } x 3 )) (2.3) 



2 

exp{-E m (p ± ) x 3 } for x 3 > 1. 



For p± = (0, 0) the long-distance behavior of these correlation functions thus defines 
electric and magnetic screening masses, which are related to the gluon polarization 
tensor, 

ml = U„(0,p 2 = -ml) . (2.4) 

Using this definition, the leading correction to the lowest order perturbation theory 
result for the electric screening mass, 



m efi (T) = J^g(T)T , (2.5) 



D Note that Gs(p±, X3) is independent of x% in Landau gauge. 



can be calculated in one-loop resummed perturbation theory. Based on the assump- 
tion that the infrared limit of the transverse gluon propagator is finite, — ^u(Po — 
0,p— > 0) = m 2 n ~ g 4 T 2 , one obtains the gauge invariant result 0, |3| 



m 2 e (T) 



m 



e.O 



Z7T 



log 



2 m P 



(2.6) 



As the magnetic mass appearing here is expected to be of 0(g 2 T) the next-to- leading 
order correction is O(glng). 

It has been pointed out already sometime ago that the leading order pertur- 
bative calculation of the electric screening mass suffers from finite cut-off effects 
I2j similarly to what has been known for bulk thermodynamic observables, e.g. 
the Stefan-Boltzmann law for an ideal gas For the Wilson action the lead- 

ing corrections are 0((aT) 2 ), i.e. 0(N~ 2 ). This is shown in Fig. |l|. For large A^ r 

1.3 

1.25 




0.07 



Figure V. The polarization tensor n o at infinite N s and various iV r . The solid line 



shows the asymptotic result (|2.7f) . 



these deviations are due to the 0(a 2 ) discretization errors introduced in the Wilson 
formulation. On a spatially infinite lattice we find for these cut-off errors 

n o(N T -> oo, iV s = oo) 2 1 f°°, 5 sinhx 



i + n; 2 — - / dxx 5 

407T 2 Jo 



2/3 g 2 T 2 ' T 40vr 2 7o sinh 4 f 
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- ^m 2 • ™ 

which is similar in magnitude to the cut-off dependence of bulk thermodynamic 



observables like the energy density ||14j| . Using an improved action these leading 
cut-off errors are eliminated and corrections only start at 0(N~ A ). In the case of 
the energy density or the pressure these actions lead to a strong reduction of cut-off 



effects in the high temperature limit WM. In the following we will see, however 



that this is not the case for the electric screening mass, i.e. the improvement of the 
ultra-violet sector does not influence the screening masses much. 

3 Determination of the Temperature Scale 

In our simulations we have used both the Wilson action (Sw) and a tree- level 
Symanzik improved action (Si), 

Sw = ijX^Trt^ixi > 

1 lxl 

Si = f^E^C/ixi-^ElV^ (3.1) 

where the sums run over all elementary lxl plaquettes £7i x i and planar 1x2 Wilson 
loops {7i x 2 |. The gauge dependent correlation functions have been calculated in 
Landau gauge, |<9 m /1 m (;e)| 2 = . Details of our gauge fixing procedure are given in 
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A comparison of simulations performed with these two actions on lattices with 
temporal extent N T = 4 will allow to judge the cut-off dependence of our results. In 
the case of the Wilson action we also have performed calculations on lattices with 
temporal extent N T = 8 in order to analyse the cut-off dependence for this action. 

In order to quantify the influence of the non-zero cut-off at finite temperature 

one should, of course, compare calculations at the same physical temperature, T = 

l/(N r a). An accurate determination of the temperature scale also is needed for 

c In the following, we will denote the couplings with f3w for the Wilson action and with /3j for 
the Symanzik improved action, respectively. 
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analyzing the dependence of the screening masses on a running coupling, g(T). We 
thus start our analysis with the determination of temperature scales for both actions. 



The problem to relate the temperature T to the coupling (3 is equivalent to the 
task of finding the dependence of the lattice spacing a on the bare coupling g 2 . We 



follow here the approach outlined in [15]. In order to take into account the violations 
of asymptotic scaling in the coupling regime of interest we use the general ansatz 



aA L 

R(g 2 ) 

48 7T 2 



R(g 2 ) ■ Kg 2 



cxp 



2 b 2 



with 
ln(& <? 2 ) - 



34/^- 
3 V167T 



(3.2) 
(3.3) 

(3.4) 



The function X(g 2 ) parameterizes the asymptotic scaling violations. For this we use 
an exponential ansatz 



Kg 2 



cxp 



^ ( d ±g 2 + ^g 4 + d zg 6 + 



(3.5) 



Using T = 1/ (N T a) we obtain from Eq. (|3.2| 



N T R(g 



Kgl) ^ 



(3.6) 



Here g 2 c is the value of the bare coupling at the critical temperature T c of the decon- 
finement phase transition at given N T . Using results for g 2 (N T ) [TjJ the function 
X(g 2 ) is obtained from a fit where T c /Al is an additional free parameter. 



Based on the Wilson action data for g 2 summarized in []16~ |, the best fit in [15 
is given by the parameterization d% = d>2 = d n>3 = 0. Their fit results are = 
5.529(63) • 10" 4 and (T C /A L ) W = 21.45(14). 

For the Symanzik improved action we have performed a similar fit, using the 
critical couplings computed in [0 for iV r > 4. Our best parameterization is given 
by d l = d n>2 = 0, and our fit results are d 2 = 5.12(18)-1(T 4 and (T C /A L ) 7 = 4.94(11). 
The fit can also be seen in Fig. ^|. 
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Figure 2: The critical temperature l/(N T R(g 2 )) vs. 4/g 2 for N T = 2,..., 8. The 
data for g~ 2 (N T ) are taken from jl7|]. The solid line is a spline interpolation of the 
data, the dashed-dotted line is obtained from a fit for N T > 4. 



Using the perturbatively calculated ratios of A-parameter, i.e. Al,i/A-l,w — 
4.13089(1) [|1|] and A m /A LjW = 19.82314 f{||, we find for the critical temperature 



a result which, within 5%, coincides with the previously determined continuum 
extrapolation for the Wilson action 

T c _ / 1.08 ± 0.01 standard Wilson action [15] , g ^ 



A^jg 1 1-03 ± 0.03 tree level improved (l,2)-action 



In the following, we will use an averaged value of T c /Aj^ = 1.06. 

We finally need to extract the temperature in units of the critical temperature 
at given N T . This is given by 

T _ R{gp ■ \{gl) 

T c R(g*) ■ \(g 2 ) • 1 • ) 



Using Eq. fl3.8|) , the fit results for A and the critical couplings from Jl6|, |7 



we 



can now relate the temperature T to the coupling (3 = 4:/g 2 . The results for the 
couplings used in our analysis are listed in Tab. [l|. The good agreement found 
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= 4 


N T 


= 8 




= 4 


Pw 


T/T c 


Pw 


T/T c 


Pi 


T/T c 


2.512 


2.004 


2.74 


2.007 


1.92 


1.984 


2.643 


3.002 


2.88 


3.031 


2.063 


3.031 


2.74 


4.013 


2.97 


3.929 


2.152 


3.923 


2.88 


6.062 


3.12 


6.016 


2.30 


5.979 


2.955 


7.527 


3.20 


7.530 


2.382 


7.528 


3.023 


9.143 


3.27 


9.151 


2.452 


9.149 


3.219 


15.88 


3.47 


15.89 


2.652 


15.88 


3.743 


66.78 


4.00 


66.71 


3.183 


66.68 


4.24 


253.5 


4.50 


253.3 


3.684 


253.2 


4.738 


953.1 


5.00 


953.9 


4.185 


954.0 


5.238 


3581 


5.50 


3578 


4.685 


3572 


5.737 


13383 


6.00 


13401 


5.186 


13393 



Table 1: Relations between the couplings and the temperatures. 

from this analysis for T c /Aj^g calculated with two different actions suggests that our 
temperature scale is of similar accuracy. 



4 Screening Masses from Gluon Correlation Func- 
tions 



In our earlier work 11 1 we have already investigated the behavior of the electric 



and magnetic screening masses in Landau gauge. Whereas in [[□]] we calculated the 
gluon propagator only at vanishing momentum, we now extend the analysis to finite 
momenta. Furthermore we are now using temperatures very much higher than in 



TT| in order to possibly get in closer contact with perturbation theory. Finally, we 
are now using in addition to the Wilson action also the tree-level Symanzik improved 
action. 

In section [| we have given the relations between the energies in the electric 
and magnetic sectors and gluonic correlation functions, Eq. ( j2.3| ). To extract the 
screening masses we use the dispersion relation between energy, screening mass and 
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momentum, which on the lattice has the form 
. -aEi . , 2 am t * . 2 a Pj 

smh = smh h k > sin , i = e, m . (4.1) 

2 2^2'' v ; 

In ( |4.1| ) we have introduced a factor k which parameterizes deviations from a free 
particle dispersion relation (k = 1) introduced by a thermal medium. 

Using T = l/(N T a) we can now compute the screening masses in units of the 
temperature, rrii/T with i = e, m. We have performed simulations using the Wilson 
action on lattices of sizes 32 3 x 4 and 32 2 x 64 x 8 and using the tree-level Symanzik 
improved action on a 32 3 x 4 lattice. At each value of the gauge coupling we per- 
formed measurements on at least 1000 configurations, see Tab. Two consecutive 



32 3 


x 4 


32 2 x 64 x 8 


32 3 


x 4 


Ay 


meas. 


Ay 


meas. 


A 


meas. 


2.512 


2000 


2.74 


1220 


1.92 


2000 


2.643 


2000 


2.88 


1000 


2.063 


2000 


2.74 


2000 


2.97 


1000 


2.152 


2000 


2.88 


2000 


3.12 


1000 


2.30 


2000 


2.955 


2000 


3.20 


1000 


2.382 


2000 


3.023 


2000 


3.27 


1440 


2.452 


2000 


3.219 


2000 


3.47 


1140 


2.652 


2000 


3.743 


2000 


4.00 


1000 


3.183 


2000 


4.24 


2000 


4.50 


1160 


3.684 


2000 


4.738 


2000 


5.00 


1000 


4.185 


2000 


5.238 


2000 


5.50 


1000 


4.685 


2000 


5.737 


2000 


6.00 


1000 


5.186 


2000 



Table 2: Number of measurement. 



configurations were separated by at least 10 update iterations, and each update 
consists of at least four overrelaxation sweeps, followed by one heatbath sweep. 

From the exponential decay of the gluon correlation functions G e , G m we extract 
the screening masses. A rather technical problem is the procedure to select a reliable 
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fit range in which G e (p±,X3) and G m (p±,X3) (see (|2.3|)) can be fitted to extract the 
energies in the electric and magnetic sectors. This is described in App. [A]. 

The results for the screening masses (from the p = measurements) and the 
energies (p ^ 0) are listed in Tabs. [3] and |] respectively. 





Wilson action, 32 


2 x 64 x 8 lattice 




Pw 


m e {T)/T 


m m (T)/T 


Pw 


m e (T)/T 


m m (T)/T 


2.74 


2.39(11) 


2.01(29) 


3.47 


1.62(4) 


0.92(7) 


2.88 


1.95(4) 


1.24(4) 


4.00 


1.62(8) 


0.66(3) 


2.97 


1.91(7) 


1.15(4) 


4.50 


1.55(5) 


0.61(2) 


3.12 


1.92(9) 


1.23(14) 


5.00 


1.41(3) 


0.52(3) 


3.20 


1.92(10) 


1.09(10) 


5.50 


1.27(5) 


0.42(2) 


3.27 


1.93(6) 


1.03(5) 


6.00 


1.26(5) 


0.37(2) 



Table 3: Electric and magnetic screening masses from G e {ki = 0) and G m (ki = 0). 



4.1 Numerical Results for Zero Momentum 

Let us first discuss the electric screening mass, extracted from the measurements 
at vanishing momentum p = 0. In Fig. |3| we show m e /T for both types of actions 
and the two different lattices we have used. One can see at once that, within errors, 
m e /T does not differ significantly for the three sets. Even the tree-level Symanzik 
improved action, which cures discretization errors of 0(a 2 ) in the action, does not 
shift the electric screening mass in any direction. This makes clear that ultra-violet 
modes do not contribute significantly to the screening mass. As a consequence, we 
have analysed all three data sets together. 

Fig. ^ shows that m e /T only depends very weakly on the temperature for small 
values of the coupling f3, corresponding to temperatures less than about 10 T c . A 
constant fit in this temperature range yields m e (T)/T = 1.938(15). This behavior 
is qualitatively similar to what we have observed in |n]]. For temperatures 1.3 T c < 
T < 16 T c we found in [|IT| m e (T)/T = 2.484(52). The difference between these 
values arises from different methods of extracting the screening masses. Whereas 
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Wilson action, 32 3 x 4 lattice 




E e (p,T)/T, extracted from 




GJki =0) 


GJki =1) 


GJki = 2) 


2.512 


2.14(11) 


2.71(18) 


2.46(5) 


2.643 


2.24(9) 


2.28(8) 


2.34(4) 


2.74 


1.94(5) 


2.13(5) 


2.33(7) 


2.88 


2.03(7) 


2.03(4) 


2.33(5) 


2.955 


1.87(4) 


1.94(4) 


2.27(5) 


3.023 


2.10(14) 


2.12(13) 


2.25(4) 


3.219 


1.80(7) 


1.93(5) 


2.03(3) 


3.743 


1.58(3) 


1.78(7) 


1.94(2) 


4.24 


1.64(8) 


1.55(3) 


2.01(4) 


4.738 


1.33(3) 


1.51(4) 


1.91(8) 


5.238 


1.19(2) 


1.40(3) 


1.83(4) 


5.737 


1.26(3) 


1.35(4) 


1.83(4) 



Symanzik action, 32 3 x 4 lattice 





E e (p,T)/T, extracted from 


Pi 


G e (k 1 = 0) 


G e (k 1 = l) 


G e {k 1 = 2) 


1.92 


2.10(6) 


2.18(5) 


2.36(4) 


2.063 


1.96(5) 


2.17(5) 


2.36(4) 


2.152 


2.08(8) 


2.04(5) 


2.23(5) 


2.30 


1.76(3) 


1.95(4) 


2.91(26) 


2.382 


2.01(11) 


2.00(4) 


2.54(9) 


2.452 


1.70(5) 


2.04(9) 


2.21(5) 


2.652 


1.72(6) 


1.53(21) 


2.05(3) 


3.183 


1.69(8) 


1.61(2) 


2.13(11) 


3.684 


1.44(8) 


1.75(6) 


1.94(2) 


4.185 


1.50(6) 


1.30(1) 


1.85(4) 


4.685 


1.19(6) 


1.45(7) 


1.74(2) 


5.186 


1.31(8) 


1.29(4) 


1.77(5) 



Table 4: Energies from the electric sector of gluon correlation functions. 
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g I i i i i i i i i l 

1 10 100 1000 10000 

T/T c 

Figure 3: Electric screening masses in units of the temperature vs. T/T c from simu- 
lations with Wilson action on 32 2 x64x8 (filled squares) and 32 3 x4 (open squares) 
lattices and with Symanzik action on a 32 3 x4 lattice (open circles). The dashed line 
is the tree- level result (|2.5|) , the dashed-dotted line is a self consistent determination 
of m e , using ( |2.6| ). The other lines are one parameter fits, using ansatz ( |4.4| ) (solid 
line, for T > 9 T c ) and ansatz ( |4~6|) (dotted line, for T > 250 T c ) respectively. 



in this paper we performed correlated fits of the gluon correlation functions over 
variable fit ranges (see App. |A]), we obtained m e and m r 



in 



111 from uncorrelated 



fits in the fixed range zT > 1. Our new method results in screening masses which 
are up to 20% smaller. Since the method accounts for possible correlations in the 
data the results should be more reliable. 



In contrast to [it] we have measured m e now also at very high temperatures (up 
to T ~ 13400 T c ; see Tab. p. From this analysis it becomes evident that m e /T 
runs with T. Since this is expected from perturbation theory it is meaningful to test 
whether perturbative predictions also work quantitatively. 

At lowest order perturbation theory for two color degrees of freedom and without 
taking dynamical quarks into account the electric mass is given by the well known 
relation ( |2.5| ). For the running coupling we use the 2-loop formula 



g- 2 {T) = T7r ^ In V + — - In 



12tt 2 A 44vr 2 



21n^ 
A 



(4.2) 
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with ijl = 2ttT being the lowest lying Matsubara frequency. Hence 



9~ (T) 



11 
12^2 



In h In 



2nT r 



A 



17 , 

+ — — o In 

447T 2 



2(ln^ + ln 27TTc 



A 



(4.3) 



We start the discussion of our data with a comparison with (|2.5| ) which is shown 
in Fig. E| as a dashed line. The numerical data for m e are lying about 60% above 



the lowest order perturbative result ( |2.5| ). However, the functional dependence of 
the electric mass on the temperature seems to be well described by m e ~ gT. 

To verify the temperature dependence of the electric mass quantitatively we 
have performed several fits of m e /T vs. ln(T/T c ) for temperatures T >9T C . In our 
one parameter fits we fix the A-parameter appearing in the temperature dependent 
running coupling to A^ and therefore use the MC-result for T c /A^, i.e. T c /A^ = 
1.06. In those cases where we parameterize the screening masses only by its leading 
g 2 dependence the effect of higher order corrections can be partially taken into 
account in a modification of the A-parameter. We, therefore, also performed two 
parameter fits with a free ratio Aflt/A^jg. 

The first fit ansatz we use is 
/ m e (T)\ 



T 



Am g\T) 



(4.4) 



The results obtained with this ansatz are summarized in Tab. [5[ They again reflect 



1-parameter fit 


2-parameter fit 


Am 

A fit/ A MS 

X 2 /dof 


1.69(2) 
1 

4.51 


Am 

A fit/ A MS 

X 2 /dof 


1.92(9) 
0.33(13) 
4.14 



Table 5: Fit results of (m e (T)/T) 2 , extracted from gluon correlation functions at 
zero momentum, using the fit ansatz ( |4.4j ). 

that the lowest order perturbative result (|2.5|) does not describe the data very well. 
The fit parameter Am is much too big compared to its theoretical value 2/3. The 
solid line shown in Fig. |3] is the result from the one parameter fit. It shows, as noted 
above, that at least the variation of m e /T with the temperature is well described by 
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ansatz ( |4.4| ). However, the temperatures we have used are apparently still too low 
to get in contact with lowest order perturbation theory. 

To test the next-to-leading order result ( |2.6|) we also determined the ratio m e /m m 
and especially the magnetic mass. We were only able to extract a reliable result for 
m m for the lattice with spatial extension N 3 = 64. On the smaller lattice the local 
screening masses m m (x3,T) do not reach a plateau (see App. Therefore the 
fits of the correlation function G m were quite poor, i.e. had a large \ 2 . As the 
electric screening masses obtained from different actions and lattice sizes do not 
show any significant difference, we expect that also the magnetic mass does not 
show a significant ultra-violet cut-off dependence. 

In Fig. [|a we show the electric and magnetic screening masses, obtained from 
the Wilson action simulation on the 32 2 x 64 x 8 lattice. Fig. |]b gives the squared 
ratio {m e /m m ) 2 . Our data strongly suggest a temperature dependence of the form 
(m e /m m ) 2 ~ g~ 2 (T), which is in agreement with the general expectation m m (T) ~ 
g 2 (T) T. We therefore performed a fit according to 



' m e {T) V 
m m (T)J 



C &t g~ 2 (T) . (4.5) 



A two parameter fit in the range T > 2 T c yields Cat = 9.16(69), Afi t /A^g = 2.42(64) 
with x 2 /dof = 0.79. Fixing the A-parameter to A^g, a one parameter fit results in 
Cfit = 7.46(27) and x 2 /dof = 1.35. In Fig. [|b we have shown the two parameter fit. 

With these results in hand we are now able to check the next-to leading order 
result for m e . The dashed-dotted line in Fig. [3] is a self consistent determination of 
m e , using (^). It lies about 20% above the lowest order prediction and therefore is 



closer to our data. However, it is still too low to describe the data well. Therefore 
we have performed additional fits of the electric mass that take into account higher 
order corrections. Based on (12.61) we use the ansatz 



+ B m g\T) . (4.6) 



' m e (T)\ 2 2 / -s/6 /m , f, 2m P 1 



log 



m m 2 



As the g A correction term leads to a temperature dependence which is too strong 
within the entire T-interval, we have restricted the fit to very high temperatures, 
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Figure 4: Electric and magnetic screening masses in units of the temperature (a) 
and squared ratio of the masses (b) vs. T/T c . Data are obtained from simulations 
on a 32 2 x 64 x 8 lattice using the Wilson action. 
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T > 250 T c . A one parameter fit at fixed T c /A m = 1.06 gives B &t = 0.744(28) with 
X 2 /dof = 4.55 (dotted line in Fig. D). 



Let us now return to the discussion of the magnetic mass. As noted above, the 
ratio m e /m m suggests a magnetic mass of the form m m (T) ~ g 2 (T)T. Therefore 
we fitted m m with the ansatz 



"im('I') _ n 2 



T 



Ait g (T) . (4.7) 



The two parameter fit of m m for T > 3 T c results in D^ t = 0.478(17) and Ag t /Aj^g = 
0.77(14) with x 2 /dof = 1.44. This is in good agreement with our result obtained 
in [ll[] for T < 20 T c . With a fixed A-parameter, T c /Aj^g = 1.06, we obtain D^ t = 
0.456(6) and x 2 /dof = 1.53. In Fig. ^a the two parameter fit is shown. The small 
deviation of the fitted curve from the measured data shows that the magnetic mass 
indeed is well described by the functional form m m (T) ~ g 2 {T) T. 

The analysis of m m /T and m e /m m is, of course, consistent with our most 
straightforward fit to the electric mass, m e {T) = ^1.69(2) g(T) T. In contrast to 
the fit based on ansatz ( |4.6|) this fit describes the data well in the entire temperature 
range above T c . This shows that the electric mass has a strong non perturbative 
character in the temperature interval we have investigated. 



4.2 Numerical Results for Non-Zero Momenta 

Let us briefly discuss the gluon correlation functions at non-zero momenta. As the 
numerical signal gets lost in statistical noise for high momenta (see ( |2.3|) and (|4.1| )) 
we only could analyse the cases k\ = 1,2, i.e. p±a = 2n/Ni, 4n/Ni. Furthermore, 
we only obtained a reliable result in the electric sector. From Eq. ( |4.1|) we have 

. 2 aEeiPl) . , 2 am e . . 2 a V\ (A ^ 

sinn = sinn h k sin . 4.8 

2 2 2 v ; 

For m e we use the result from the measurement at zero momentum. In the limit 
T — > oo one expects to find a free particle dispersion relation, i.e. k — > 1. In Fig. [| 
we have plotted k vs. T/T c . Obviously we do not have sufficient statistics to uncover 
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Figure 5: k vs. T/T c for k\ = 1, 2. Some data points have been displaced horizontally 
for better viewing. 

a temperature dependence of k. Therefore we only quote a value averaged over the 
temperature interval T >9T C . We find k = 0.37(10) for k\ — 1 and k = 0.65(3) for 
k± = 2. This suggest a quite significant modification of the free particle dispersion 
relation at low momenta. 



5 Polyakov Loop Correlation Functions 

It is well known that for temperatures above the critical temperature T c , the confine- 
ment potential between a quark and an anti-quark is replaced by the color averaged 
potential p0| , which, in lowest order perturbation theory, is of the form 

V m (R,T)~^e- 2m °W R for T > T c . (5.1) 

Here m e is the electric (or Debye) screening mass. As V av decreases very fast, the 
numerical signal gets lost in statistical noise in the long distance regime. On the 
other hand, ( |5.1| ) is only valid at large distances. This situation is improved for the 
color singlet potential, which to leading order perturbation theory is controlled by 
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1-gluon exchange and therefore takes on the form 



N 2 _ i -m e (T)R 

The color singlet potential, however, is gauge dependent and one therefore again has 
to fix a gauge for its calculation. As mentioned earlier we have chosen the Landau 
gauge. 

On the lattice one can extract both potentials by measuring Polyakov loop cor- 
relation functions pOi . 



p -VMR,T)/T _ (TrL(fi)TrLt(O)) 

-Vi(R,T)/T _ 

o (Tr(L(i?)Lt(0))) 

e - (|L|)2 . ^.4j 



5.1 Improvement of the Singlet Potential 

In this section we want to discuss briefly the improvement of rotational symmetry for 
the color singlet potential due to the use of an improved action. We have measured 
V\/T both along an axis, labeled with (1,0,0), and along three different off-axis 
directions, (1,1,0), (1,1,1), and (2,1,0). To make the results from simulations 
with unimproved and improved action comparable, one has to choose couplings that 
both correspond to the same temperature. As an example, we use in the following 
Pw — 3.219 and /3j = 2.652. As listed in Tab. [H, both couplings correspond to 
T ~ 15.88 T c . 

Motivated by Eq. ( |5.2j ) and taking into account the periodic boundary conditions, 
we have performed a correlated fit of the (1, 0, 0) data in the interval^ R 6 [7, 12], 
using the fit function 

M«)=4(— (5.5) 

with A^3 = 32. The fit results for both actions are listed in Tab. |6|. As one can see 
d See A PP . 0. 
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fit parameters 


p w = 3.219 


Pi = 2.652 


A &t 


-1.49(16) 


-1.64(13) 


met 


0.435(16) 


0.416(11) 


goodness 


0.380 


0.452 


X 2 /dof 


1.049 


0.918 


normalized x 2 deviation from the (1,0,0) fit 


(1,1,0) 


2.019 


1.457 


(1,1,1) 


1.934 


0.090 


(2,1,0) 


1.316 


0.243 



Table 6: Results from the fits of V x /T at (3 W = 3.219 and /3j = 2.652 (T ~ 15.88 T c ). 



from the upper part of the table, the fit itself is better for the improved data than 
for the Wilson data, i.e. the errors on the fit parameters are smaller, the goodness 
is higher and finally the squared error from the correlated fit (x 2 /dof) is smaller. 
The lower part of Tab. ^| shows the \ 2 deviation of the off-axis data points from the 
(1, 0, 0) fit curves. For this comparison we used data in the interval 7 < R < 12 and 
divided by the number of points taken into consideration. For all measured off-axis 
directions these data show that the violation of rotational symmetry is lowered by 
going from the Wilson to the tree-level Symanzik improved action. This behavior 
becomes also clear from Fig. |6|, which shows the potential V\ (normalized by the fit 
function ( p.5| ) with the parameters given in Tab. []) vs. distance R . 



5.2 The Electric Screening Mass from the Singlet Potential 

We used on-axis point-to-point as well as plane-plane Polyakov loop correlation 
functions to extract the electric screening mass. In the former case, we used Eqs. 
( |5.2|) and ( |5.4|) and performed a correlated fit of the numerical data, using ( |5.5| ) and 
the fit criterium described in App. We did this both for the measurement along 



the (1,0,0) axis and for the three different off-axis directions mentioned in Sec. [54 



In the second case, the expression for the Polyakov loop L(R) in Eq. ( |5.4|) is 
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Figure 6: Singlet potential Vi(R), normalized by the correlated fit of the (1,0,0)- 
data in the interval R G [7, 12]. The data have been calculated on a lattice of size 
32 3 x 4 using the Wilson action at (3w = 3.219 (a) and the Symanzik action at 
Pi = 2.652 (b). Both couplings correspond to a temperature of T ~ 15.88 T c . The 
different symbols refer to the (xi, X2, xs)— directions along which the measurements 
have been performed. 
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Figure 7: Electric screening masses, obtained from Vi )Sum . The data points refer to 
the same lattice sizes and actions as in Fig. |3|. The solid line is a one parameter fit 
for T > 9 T c , using ansatz 



replaced by L(x 3 ) = J2 Xl ,x 2 L(xi,X2,x 3 ). Then (|5.4j ) and (|5.2|) transform into 



e -V llBwa {x 3 ,T)/T _ 2 



(Tr(L(x 3 )Li(0))) 
(\L\) 2 



(5.6) 



and 



yX 3 



T) ~ e _me(T)xa for T > T c 



(5.7) 



Whereas we have calculated V\ on lattices of size 32 3 x 4 and 32 2 x 64 x 8 and for 
both actions, we have calculated Vi )Sum only on the larger lattice, using the Wilson 
action. The results for the electric screening mass are listed in Tab. [7] and || 

Similar to the electric mass extracted from gluon correlation functions, the results 
we have now obtained with different actions and on lattices of varying size again 
do not differ significantly. Therefore we have also here analysed all three datasets 
together. The screening masses, extracted from Vi jSum , are shown in Fig. [F[ 

As expected from Sec. [|, m e /T, extracted now from Polyakov loop correlation 
functions, also depends only weakly on the temperature for temperatures less than 
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Wilson action, 32 3 x 4 lattice 







m e (T)/T, extracted from 




Pw 


\A 

v l,sum 


r 1,(^1, U,UJ 


l/i /i -1 ^^ 
* 1,(1,1,1); 


K 1,(1, Ij 1 ; 


•'1,(2,1,0) 


9^19 
Z.OIZ 


9 rt'\( r >\ 

z.uo(z j 


2 28(9) 


2.22(4) 


2.13(4) 


Z.oo(l 1 J 




z.ou^y j 


2.15(6) 


2.12(5) 


2.07(4) 


9 9^0"! 

z.zo(y j 


9 7/1 
Z. 1 4t 


z. io(y ) 


2 09(5) 


2.04(5) 


2.04(4) 


9 1 7^ 
Z. l^ 1 j 


9 88 
Z.oo 


9 D/1 ffi"! 
Z.U^O ) 




1 Q1 (^ 




1 QR^7^ 

i.yo( i ) 


once 

2.955 


1.94(7) 


2.11(8) 


1.97(5) 


1.94(6) 


2.06(5) 


3.023 


2.16(7) 


1.84(5) 


2.01(4) 


1.96(9) 


1.98(6) 


3.219 


1.88(9) 


1.74(6) 


1.80(4) 


1.75(5) 


1.83(4) 


3.743 


1.85(15) 


1.50(2) 


1.50(1) 


1.43(2) 


1.61(3) 


4.24 


1.63(7) 


1.55(9) 


1.53(7) 


1.45(3) 


1.44(2) 


4.738 


1.33(4) 


1.75(10) 


1.32(5) 


1.32(3) 


1.29(2) 


5.238 


1.30(7) 


1.29(4) 


1.25(4) 


1.20(2) 


1.19(2) 


5.737 


1.34(5) 


1.29(3) 


1.24(4) 


1.17(1) 


1.15(2) 



Symanzik action, 32 3 x 4 lattice 







m e (T)/T, extracted from 




Pi 


^l,sum 


^1,(1,0,0) 


Vi, (1,1,0) 


^1,(1,1,1) 


^1,(2,1,0) 


1.92 


2.26(9) 


2.04(4) 


2.16(5) 


2.12(5) 


2.06(2) 


2.063 


1.97(3) 


1.98(3) 


2.01(3) 


2.03(4) 


2.12(7) 


2.152 


2.07(5) 


2.16(8) 


2.02(7) 


2.04(4) 


2.01(6) 


2.30 


1.93(5) 


1.83(2) 


1.79(3) 


1.70(3) 


1.84(4) 


2.382 


1.82(4) 


1.83(3) 


1.72(3) 


1.70(2) 


1.90(5) 


2.452 


1.78(6) 


2.14(12) 


1.98(14) 


1.76(7) 


1.79(5) 


2.652 


1.69(5) 


1.66(4) 


1.57(2) 


1.54(2) 


1.67(3) 


3.183 


1.65(6) 


1.58(6) 


1.62(5) 


1.51(3) 


1.48(2) 


3.684 


1.38(4) 


1.71(12) 


1.42(3) 


1.36(3) 


1.33(2) 


4.185 


1.62(7) 


1.66(15) 


1.34(3) 


1.19(2) 


1.21(66) 


4.685 


1.46(12) 


1.22(3) 


1.19(3) 


1.10(2) 


1.09(1) 


5.186 


1.16(2) 


1.22(5) 


1.14(3) 


1.05(2) 


1.03(1) 



Table 7: Electric screening masses from Polyakov loop correlation functions. 
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Wilson action, 32 2 x 64 x 


8 lattice 


Pw 


m e (T)/T 


Pw 


m e {T)/T 


Pw 


m e (T)/T 


2.74 


2.46(16) 


3.20 


1.76(4) 


4.50 


1.55(10) 


2.88 


1.89(4) 


3.27 


1.93(5) 


5.00 


1.36(3) 


2.97 


1.80(31) 


3.47 


1.65(3) 


5.50 


1.18(8) 


3.12 


1.92(8) 


4.00 


1.61(4) 


6.00 


1.16(31) 



Table 8: Electric screening masses from Vi jSum . 

about 9T C . For example, fitting Vi )Sum in this temperature range with a constant, 
we obtain m e (T)/T = 2.010(13). 

According to Sec. |] we performed, for temperatures T > 9 T c , one parameter 
fits (i.e. Afi t /Aj^g = 1), using the ansatz Q4.4Q . The results from the lattice of size 
32 3 x 4 are summarized in Tab. |9|. On the 32 2 x 64 x 8 we obtain from Vi iSum a fit 



Wilson action, 32 3 x 4 lattice 




Fit 


s of (m e (j 

^1,(1,0,0) 


0/T) 2 , e> 

^1,(1,1,0) 


:tracted fr 
^1,(1,1,1) 


om 

Vl, (2,1,0) 


Ait 

X 2 /dof 


1.97(6) 
2.49 


1.70(3) 
10.6 


1.60(2) 
6.64 


1.53(2) 
7.29 


1.61(2) 
4.18 



Symanzik action, 32 3 x 4 lattice 




Fit 


s of (m e (j 
Vi, (i,o,o) 


n/Tf, e> 

^1,(1,1,0) 


:tracted fr 


om 

Vl, (2,1,0) 


X 2 /dof 


1.67(4) 
7.10 


1.46(2) 
15.0 


1.46(2) 
10.8 


1.35(2) 
4.51 


1.36(2) 
5.72 



Table 9: Fit results of (m e (T)/T) 2 , extracted from Polyakov loop correlation func- 
tions, using the fit ansatz ( |4.4| ). 

value Afi t = 1.72(4) with x 2 /dof = 4.60. 

In general we find that the results extracted from Vi iSum are in good agreement 
with the zero momentum results from the gluon correlation functions. To make this 
clear also quantitatively we have analyzed all three datasets for Vi )Sum together, as 
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in the case of the gluon correlation functions. The one parameter fit for T > 9T C 
yields Am = 1.71(2) with x 2 /dof = 5.80. This can be compared with the result from 
Tab. §, v4 fit = 1.69(2) with x 2 /dof = 4.51. We therefore conclude that the electric 
screening mass is well described by m e {T) = ^1.70(2) g(T) T in the temperature 
range T < 14000 T c . 

6 Summary and Conclusions 

We have studied Polyakov loop and gluon correlation functions in the high tempera- 
ture deconfined phase of SU(2) lattice gauge theory in a wide range of temperatures, 
using both the standard Wilson action and a tree-level Symanzik improved action. 
We have calculated chromo-electric and -magnetic screening masses and have deter- 
mined their dependence on the temperature. 

The temperature dependence found for the magnetic mass is in accordance with 
the expected g 2 (T)-dependence. We find m m {T) = 0.456(6) g 2 (T) T. The behav- 
ior of the electric mass shows, however, that this does not at all mean that the 
screening masses can be described according to perturbative expections. Although 
the temperature dependence of m e is consistent with a logarithmic dependence, 
m e ~ gT, our data do not agree with lowest order perturbation theory. Only little 
improvement is achieved by using next-to-leading order results from resummed PT. 
From an analysis of the gluon propagator as well as the color singlet potential we 
find m e {T) = ^1.70(2) g(T) T. This result shows that the screening mechanism is 
highly non-perturbative even for temperatures as large as 14000 T c . This observation 
is in accordance with studies of screening in dimensionally reduced 3d-QCD || ||. 

Our simulation of the gluon correlation functions at finite momenta still suffer 
from insufficient statistics. We find a modification of the energy momentum disper- 
sion relation of a free particle, but we are not yet able to quantify its temperature 
dependence. 

The improvement of the action does not show, within statistical errors, any 
significant modification of the behavior of the screening masses, although we can 
show that the violation of the rotational symmetry of the singlet potential, which 
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also was used to extract m e , is weakened. 
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Appendix 

A Determination of Screening Masses from Cor- 
relation Functions 



To shorten the discussion, we assume in this section the case of vanishing momen- 



and (|4.1|) that the electric and magnetic screening masses are related to the long- 
distance behavior of the gauge field correlation functions G e (x 3 ) = G e (p± = 0,x 3 ) 
and G m (x 3 ) = G m (p± = 0,x 3 ) via 



As there is generally no unique rule how to select the interval in which G e (x 3 ) and 
G m {x 3 ) should be fitted, we have constructed a criterium how to find a suitable fit 
interval for a given set of numerical data. 

A fit is in general characterized by several properties, namely x 2 , goodness (Q), 
degrees of freedom (g) and relative errors of the fit parameters. As x 2 enters directly 
into the calculation of Q, it is sufficient to consider only the last three quantities. 
It is desirable to find a fit with large Q, large g and small relative error on the 
fit parameter we are interested in, i.e. on rrii/T. As we want to weigh these three 
quantities, we are looking for a fit interval with 



turn. Taking periodic boundary conditions into account, it follows from Eqs. ( |2.8| ) 




for x 3 ^> 1 and i = e,m. 



(Al) 




(A.2) 
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We have chosen the coefficients to be a = 9, (3 = 1 and 7 = 3. As we put the largest 
weight on Q, it sometimes happens that only a very small fit interval is selected 
by this condition. To avoid this problem, we add the extra condition g > 3. For a 
two parameter fit this is equivalent to demand that the fit interval should contain 
at least 5 points. 

Finally we require that the fit only considers points G e (x 3 ) or G m {x 3 ) with errors 
less than 50 % of their value. This is to reject points that are dominated by statistical 
noise. 

To see how our fit criterium works we show in Fig. |8| G e (x 3 ) vs. x% for one arbi- 
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Figure 8: The electric correlation function G P { X3) clS cL function of £3 calculated on 
a lattice of size 32 2 x 64 x 8 using the Wilson action at flw = 3.47. The squares 
show points with an error less than 50% of the value, whereas the circles describe 
points with a bigger error. The filled points represent the fit interval, found by our 
fit criterium. The solid line is the correlated fit. 

trary coupling, f3 w = 3.47. The solid line is the correlated fit, found automatically 
by the fit criterium. To demonstrate the quality of the fit criterium, we have also 
studied local screening masses m^Xs). They are defined by the relation (a = 1) 



Gi(x 3 ) cosh(mi(x 3 ) 



X3 2 3 



Gi(x 3 + 1) cos h (mi(x 3 ) (x 3 + 1 - ^ 



e, m 



(A.3) 
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If £3 becomes large enough, rrii{x^) must reach a plateau. On the other hand, if 
£3 becomes too large, the local masses have big statistical errors and do not carry- 
valuable information. 



In Fig. |9| we show the local electric screening masses in units of the temperature, 




Figure 9: Local electric screening masses, extracted from the electric correlation 
function G e (xs) shown in Fig. || The horizontal lines are the lower and upper 
bounds for m e (T)/T, given by the correlated fit shown in Fig. ||. 

m e (x 3 ,T)/T, extracted from G e {x^) shown in Fig. |^. The horizontal lines are the 
lower and upper bounds for m e (T)/T, given by the correlated fit shown in Fig. 
Obviously, the fit criterium works fine and finds a reliable fit interval. 
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